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It has long been hypothesized that futile cycles in cellular metabolism are involved in the regulation 
of biochemical pathways. Following the work of Newsholme and Crabtree, we develop a quantitative 
theory for this idea based on open-system thermodynamics and metabolic control analysis. It is 
shown that the stoichiometric sensitivity of an intermediary metabolite concentration with respect 
to changes in steady-state flux is governed by the effective equilibrium constant of the intermediate 
formation, and the equilibrium can be regulated by a futile cycle. The direction of the shift in 
the effective equilibrium constant depends on the direction of operation of the futile cycle. High 
stoichiometric sensitivity corresponds to ultrasensitivity of an intermediate concentration to net flow 
through a pathway; low stoichiometric sensitivity corresponds to super-robustness of concentration 
with respect to changes in flux. Both cases potentially play important roles in metabolic regulation. 
Futile cycles actively shift the effective equilibrium by expending energy; the magnitude of changes in 
effective equilibria and sensitivities is a function of the amount of energy used by a futile cycle. This 
proposed mechanism for control by futile cycles works remarkably similarly to kinetic proofreading 
in biosynthesis. The sensitivity of the system is also intimately related to the rate of concentration 
fluctuations of intermediate metabolites. The possibly different roles of the two major mechanisms 
for cellular biochemical regulation, namely reversible chemical modifications via futile cycles and 
shifting equilibrium by macromolecular binding, are discussed. 

PACS numbers: 



I. INTRODUCTION 

Metabolic fluxes and intermediary metabolite con- 
centrations are the essential currency of intracellular 
metabolism. The fluxes — or rates of turnover of the var- 
ious reactions — determine the rate of production of re- 
quired end products of a metabolic system, i.e., ATP 
synthesis by oxidative phosphorylation and lactate pro- 
duction via glycolysis. The intermediary metabolite con- 
centrations are often key regulators in metabolic con- 
trol. For example, fructose- 1,6-bisphosphate is an al- 
losteric regulator of pyruvate kinase . One of the pri- 
mary examples of futile cycles in metabolism is fructose- 
6-phosphate phosphorylation dephosphorylation^n ex- 
tensively studied metabolic regulatory module 
In the present work, we study the quantitative relation- 
ship between a metabolic flux and the concentrations of 
its intermediary metabolites. In particular, we study the 
sensitivity of the intermediate concentration with respec- 
tive to changes in the flux. Depending on the biological 
context, the sensitivity may be high for some metabolites 
(sensitive) while for others it may be low (robust) . We 
show that the sensitivity is intimately related to the en- 
ergetics, i.e. effective equilibrium constant, between the 
intermediate and the upstream substrate(s). Hence con- 



* Electronic address: "qian Oamath. Washington. edu| 



trol of sensitivity has a thermodynamic interpretation. 

The chemical equilibrium constant is an intrinsic prop- 
erty for a given chemical reaction and cannot be altered 
by enzyme activity. However, the effective equilibrium 
constant operating in a reaction in a network can be 
shifted via a futile cycle, as we describe below. This 
same idea is behind the concept of kinetic proofread- 
ing, in which a kinetic cycle involving GTP hydrolysis 
increases the effective binding affinity between a codon 
and its tRNA, improving the accuracy of protein biosyn- 
thesis 0, S 13 • At a deeper level, the kinetic proofread 
shares a same principle of nonequilibrium physics as the 
nuclear Overhauser effect in magnetic resonance and cat- 
alytic wheel 0- 

The idea that futile cycles contribute to the regulation 
of metabolic functions has a long history 0,l3,llO, HJ- In 
particular, Newsholme and Crabtree suggested that fu- 
tile cycling may be important for effective regulation at 
low metabolic fluxes, allowing a pathway to be controlled 
with much smaller excursions in the concentrations of 
the allosteric regulators. However, there has not been a 
systematic mathematical theory substantiating this idea. 
Because of the lack of quantitative predictions, experi- 
mental verification of this important idea has been diffi- 
cult. The present analysis provides strong support and 
a mechanistic quantification of the hypothesis of New- 
sholme & Crabtree on the biochemical function of futile 
cycles in metabolic systems [l^, E3 • The analysis focuses 
on the quantitative relation between the energy expendi- 
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ture, thermodynamics, and biochemical control. 



II. METABOLIC FLUX, INTERMEDIATE 
CONCENTRATION, AND STOICHIOMETRIC 
SENSITIVITY 

A. Sensitivity of Intermediate Concentration to 
Flux through a Metabolic Pathway 

Consider a simple linear metabolic pathway with 
species 0, 1, and 2: 



k_i k_ 



(1) 



in which the concentration of the substrate, cq, is fixed. 
The constants ki are either the mass-action rate con- 
stants for first-order chemical kinetics or the effective rate 
constants acting around a steady-state. (See Appendix 
A for a detailed analysis.) Note that in a metabolic net- 
work, there are different ways of maintaining the input 
substrate concentration. Two extreme cases are concen- 
tration clamping and constant flux injection Q|. While 
concentration clamping can be achieved experimentally 
with a large pool or buffer for the substrate, flux injection 
is closely related to the "flux-generating step" suggested 
in [l5| . Both concentration clamping and flux injection 
provide thermochemical driving forces for open biochem- 
ical systems, analogous to batteries in electrical circuits. 

Recall that there are two types of ideal batteries, those 
that provides constant voltage (voltage sources with zero 
internal resistance) and those that provides constant cur- 
rent (current sources with zero internal conductance). 
A real battery of course has a finite internal resistance 
and conductance. Metabolic concentration and flux play 
equally important roles in the steady-state of a biochem- 
ical system. They should be treated on equal footing 
in a complete theory of metabolic dynamics. Under an 
ideal setting, one can control the concentration(s) and 
let the fluxes change in response; Similarly, one can con- 
trol the flux(s) and let the concentrations change in re- 
sponse. Controlling fluxes can, but not necessarily, be 
accomplished by changing enzyme activities. 

Significant controversies exist in the literature on sys- 
tems analysis of metabolic networks due to differences in 
implicit assumptions on how a system is sustained in a 
nonequilibrium steady-state. In other words, how a sys- 
tem's steady-state responses to perturbations depends on 
how the steady-state is maintained. See T^l for an inter- 
esting case study. It is worth pointing out that almost all 
the existing work on metabolic control analysis (MCA) 
implicitly assumes concentration clamping [17i| . In fact, 
we recently discovered that for systems driven by flux in- 
jections, the summation of flux control coefficients equals 
zero rather than unity (manuscript in preparation) . Nat- 
urally, the realistic situation in a cell is likely to be mixed. 

In standard MCA, one defines a flux control coefficient 
as the sensitivity of a flux J to an enzyme concentration 



[TM IT^ l20l|. This work, however, focuses on a different 
aspect of network control: the steady-state sensitivity of 
the concentrations of the intermediates, ci and C2, to the 
flux J, at fixed enzyme activities. This stoichiometric 
sensitivity is also different from the elasticity coefficient, 
which is a property of a single enzyme — an elasticity co- 
efficient (also called local T7^ or intrinsic TFl, or im- 
mediate ,21,| control coefficient) is determined from the 
rate law of a single enzymatic reaction in isolation. Our 
stoichiometric sensitivity is related to the co-response co- 
efficients introduced by Rohwer and coworkers [l^, 113 ■ 
The co-response coefficient emphasizes the concomitantly 
changes in a flux and a concentration in response to a per- 
turbation in a given enzyme. While changing enzyme is 
one of the possible means to perturb kinetics, there are 
other means to perturb a flux. Hence the stoichiomet- 
ric sensitivity is a more general kinetic concept than the 
co-response coefficient in MCA. In fact, the present anal- 
ysis focuses on sensitivity of network concentrations to 
steady-state network fluxes. Even though both metabo- 
lites and enzymes are chemical species in a biochemi- 
cal reaction system, their roles in metabolic kinetics and 
thermodynamics are very different. A change in the en- 
zyme amounts to a change in both forward and backward 
effective rate constants for the catalyzed reaction without 
altering their ratio (Haldane's equation). 

With given cq and J as input and output in a steady- 
state, the intermediate concentrations in the reaction of 
Eq. n are 



Cl 



C2 



A. 

fc_i 



Co 



J 



kik2 

fc-lfc_2 



Co 



k2 + k^i 

fc-lfc-2 



(2) 



(3) 



(See Appendix A.) The stoichiometric sensitivity coeffi- 
cients are defined as 



'71 



V2 




J 



(fca + fc-i)J 

fc_ifc_2C2 



(4) 



(5) 



While the steady-state concentations do not necessarily 
increase or decrease in the order from input to output 
along the pathway, the sensitivity increases for interme- 
diates as one moves from "upstream" to "downstream" 
(??2 > rji since fc2Ci — k-2C2 > 0). This observation gen- 
erally holds true. For a sequence of reaction of arbitrary 
length 



^ n 



(6) 



the stoichiometric sensitivity for the ith intermediate is 

I -\- fcj _|_ _|_ ki---k2 

= ^-"-^» , ^zllzi)l±^ J. (7) 
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This expression is related to the well-known exit proba- 
bility and mean first passage time out of state i |2J| , as 
is shown in . The Gibbs free energies of the metabo- 
lites, Gi = G° + RTlnci, decreases along the pathway 
following the direction of the flux. As we shall see, this 
correlation between the sensitivity and energetics is not 
just a coincidence. 

The sensivity rji can also be expressed as J/ J^, where 
J_ is the backward flux from state 1 to state 0, and the 
net flux J is equal to the foward flux minus the backward 
flux J ~ — J_. I I = J_ / J is an elegant result 
that was known to Newsholm and Crabtree The 
flux ratio is a measure of whether a reaction operates near 
equilibrium (J/ J- <C 1) or far from equilibrium (J/ J- 3> 
1). The expression rji = J/ J_ applies to an intermediate 
at an arbitrary position in the reaction sequence when 
all the upstream reactions are in rapid equilibrium. For 
example, if 3> k2, then the reaction between state 
and state 1 is maintained near equilibrium, and 772 = 
J/ J- for the reaction between state 1 and state 2. 

Quantities such as gj" ^ have been discussed in the 
context of allosteric regulation, by a metabolite c, of 
some enzyme which in turn regulates the J Brand 
and his coworkers have extensively used an empirical, 
"top-down" elasticity analysis in assessing the fractional 
changes in the fluxes of metabolic reactions in response 
to a change in the concentration of an effector [26l| . 
These studies, which emphasize the interactions between 
metabolites and enzymes, are different from the above 
direct sensitivity of J with respect to an intermediate 
concentration due solely to the nature of stoichiometric 
networks. 

In terms of the stoichiometric sensitivity, we now ask 
a typical engineering question. How can one reduce (or 
increase) the sensitivity of ci with respect to J, with 
given substrate source cq, flux J, intermediary metabolite 
concentration ci? We assume here that ci is a regulator 
and/or control agent for some other parts of the cell, 
the specific nature of which is not of our concern. To 
change iji according to Eq. ^ while maintaining ci at 
its fixed value, the effective fci and must be changed 
simultaneously so that 

which leads to 

^ _ / Sk_i \ cifc_i 

fci V fc_i / Cofci 

However, by altering enzyme activity, the effective fci 
and fc_i can be changed only according to Ski/ki — 
dk-i/k-i. Thus the sensitivity around a set-point flux 
and intermediate concentration cannot be effectively ad- 
justed by increasing or reducing the enzyme activity for 
the reaction cq ^ ci. To adjust fci and fc_i indepen- 
dently of one another it is necessary to change the ef- 
fective AG° = i?T ln(fc_i/fci) for reaction 0^1. (Do- 
ing this is known as impact control in \2T\). From Eqs. 



(0) and it can be shown that J = fcico — fc_ici = 
J_(fciCo/fc-iCi — I), hence 

Cl 

Eq. ^1 is our flrst key result. It shows that the sensi- 
tivity of an intermediary metabolite concentration with 
respect to the steady-state flux is governed by the AG°, 
the Gibbs free energy of the intermediate formation. Fig. 
1 shows how 771 changes as a function of the AG" . Qual- 
itatively, this can be understood without the mathemat- 
ics: There are two extreme cases of a steady state of the 
pathway in Case I (far from equilibrium) is that J = 
fciCQ — fc-ici fc_ici. In this case AG°/i?r ^ ln(co/ci) 
and 771 3> 1. Case II (near equilibrium) is that the fc_ici 
3> J. In the case where the reaction ^ 1 is maintained 
near equilibrium, there is little sensitivity of ci to flux J. 
High sensitivity occurs when fci fc_i; all other things 
being equal, the smaller the value of fc_i, the smaller the 
value of AG°, and the greater the sensitivity. Note that 
771 goes to zero when cq and ci are in equilibrium; 771 
computed by Eq. ^1 remains positive for all cases for 
which the reaction proceeds in the positive direction. 



60 




FIG. 1: Relation between Gibbs free energy, AG° and the 
stoichiometric sensitivity 771. For reaction 1 with 
co/ci — 8. Smaller the AG°, greater the sensitivity (771) of 
Cl in response to output flux (Eq. llOt . Everything else being 
equal, greater AG° means greater backward flux from ci to 
Co. When it is significantly greater than the net flux J, the 
level of Cl will be insensitive to changes in J. 



B. Impact of Futile Cycles on Stoichiometric 
Sensitivity 

For isolated chemical reactions, modifying AG° can be 
accomplished only through modifying the solvent condi- 
tions. Such a mechanism is clearly not of primary impor- 
tance for cellular biochemical reactions. Structural mod- 
ifications change the nature of the chemical reactions. 
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leading to different values of AG°. This mechanism rep- 
resents a possible biological strategy from an evolution- 
ary standpoint, but it is not useful as a mechanism for 
dynamic regulation of the sensitivity in a cell. 



i?rin 



ki 



is the chemical driving force in the futile 



(«) 




FIG. 2: Shown in (a), a biochemical reaction between specie 
and 1 in isolated system reaches its equilibrium with concen- 
trations Icq = ki/k-\ = e~^^ IRT ^ Enzyme can change 
the rate constants, but not the free energy difference AG°. 
However, if this reaction is coupled to other reactions in an 
open biochemical network as shown in (6), a futile cycle is 
able to shift the population ratio ci/co to be greater (or less) 
than the equilibrium value fei/fc_i. In (&), the additional re- 
actions involve species D and E. There is now a futile cycle 
involving species and 1. The equilibrium between D and 
E is c^jc^ ~ Ci''fc3/(cQ''fc_3) — k\k-i,l(k-\k-2). If the con- 
centrations of D and E are not at their equilibrium, then 
XnicEkxkz I {cDk-xk-z)) = AGde 0, which is the active en- 
ergy source (e.g., nucleotide hydrolysis) that pumps the futile 
cycle. In a steady state this energy is dissipated as heat. The 
same mechanism is behind the nuclear Overhauser effect in 
magnetic resonance, kinetic proofreading in biosynethesis 0, 
and catalytic wheel 0. 

Since the free energy of formation of an intermediary 
metabolite (at a given temperature, pressure, and sol- 
vent condition) cannot be altered, is there a solution to 
reducing (increasing) the sensitivity in a reaction net- 
work? One of the possible mechanisms for increasing or 
reducing 771 is a coupling between the reaction and a fu- 
tile cycle. Fig. 2 shows a futile cycle attached to the 
reaction 0^1, where if '^"^"^^^'^ ^ 1, then the appar- 
ent free energy difference between species and 1, AG" 
is 



^AG°/RT 



in which a 



ki + fc_3 



1 



^AGde/RT 



1 



k^CE and fc_ 



(11) 

k-3CD 



cycle in Fig. 2. By apparent free energy, we mean that 
we would treat the nonequilibrium steady-state concen- 
tration ratio between and 1 as they were in an equi- 
librium: ci/co = e~^'^ /RT_ ^] the sensitivity 
771 is expressed as a function of AG". With the reaction 
coupled to the futile cycle illustrated in Fig. 2, AG° 
AG°, and rji is expressed as a function of AG° 



coe 



-AC/RT 



Cl 



Eq. El can be expressed 



1 



(12) 



(13) 



where 77° is equal to ryi in the absence of the futile cycle 
(AGdb = 0). This is our second key result. It shows that 
the change in sensitivity rji can be controlled through the 
amount of available energy, i.e. the driving force AGde 
(also see Appendix B). The amount of energy consumed 
by each turn of the futile cycle is \AGde\- Fig- 3 shows 
how the sensitivity rji varies as functions of the energy 
expenditure AGjje for different tr, the relative rates of 
the two steps in the loop. 
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FIG. 3: Stoichiometric sensitivity (771) as a function of the 
driving force in the futile cycle (AGde) given in Eq. (1131 . 
with different values for the parameter a: (a) a — 10, (b) a = 
1, and (c) a — 0.01. rj" is the corresponding rji when AGde ~ 
0, i.e., without the regulation from the futile cycle. Note from 
Fig. 2, when the AGde < the flux in the futile cycle goes 
clockwise and when AGde > 0, it is counter clockwise; They 
correspond to 771 < 77° and 771 > rjf, respectively. 

A third key result follows from Eq. ^1 which allows 
us to diagnose the putative role of a futile cycle in a bi- 
ological system. When AGde < 0, then the futile cycle 
of Fig. 2 is driven in the clockwise direction, moving the 
reaction (at fixed cq and ci) away from equilibrium, and 
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FIG. 4; Futile cycle (a) versus macromolecular binding (b), 
two mechanisms for modifying the equilibrium concentra- 
tion ratio between species and 1. (a) The futile cycle 
scheme requires active energy expenditure AG but only small 
amount of proteins as the enzyme for the additional reac- 
tion. fc±3 are pseudo-first order rate constants: fca = k^CE 
and fc_3 = k-iCD in Fig. 2. When AG = 0, the population 
ratio becomes Keq- {h) The macromolecular binding scheme 
requires no energy cost, but it requires a large amount of pro- 
tein being synthesized in advance. When the concentration 
of if the protein (P) Cp = 0, one has AG* = — oo and the 
population ratio becomes K^q- 



increasing the sensitivity of ci to changes in flux. If the 
futile cycle is thermodynamically driven in the counter- 
clockwise direction, i.e., when /S.Gde > 0, then sensitiv- 
ity of ci is reduced, i.e., the intermediate concentration 
is made robust to changes in flux. In both cases there is 
an energy expenditure and heat dissipation. The "high 
grade" chemical energy is transformed into "low grade" 
heat, but is not wasted from the standpoint of informa- 
tion regulation. 



C. Sensitivity of Intermediate Concentration to 
Input Concentration 

To investigate the potential impact of a futile cycle on 
the sensitivity of intermediate concentration to the input 
concentration, we define Ci as the sensitivity of interme- 
diate concentration ci to changes in input concentration 
Co, at a given steady state flux J, 



Ci 



31nci 
d In Co 



(14) 



From the equation ci = e ^g°/rt^^ _ J/(fciCo))co, it is 
straightforward to show that 



(15) 



when there is no futile cycle acting on the reaction. For 
the case when the futile cycle of Fig. 2 is present. 



cij Vl + cre 



AGl 



(16) 



Equation El shows that the futile cycle has qualitatively 
the same impact on intermediate concentration sensitiv- 
ity to input concentration as on intermediate concentra- 
tion sensitivity to flux. 



III. OPTIMAL INTERMEDIARY METABOLITE 
SENSITIVITY AND ROBUSTNESS 

The previous discussion was based on the example of 
Eq. We now present a more general theory for sto- 
ichiometric sensitivity in biochemical networks. For a 
given reaction r ^ p in a biochemical network, its kinet- 
ics, according to the law of mass actions, and thermo- 
dynamics are determined by its forward and backward 
fluxes J+ and J_ 



J=J+-J_, AG = i?Tln 



J-_ 



(17) 



If we perturb concentrations and flux around a steady- 
state we have 



6 J = SJ+ - SJ-, 

dJ- SJ+ 
J- J+ 



SAG 



RT 



(18) 



Solving (5J+ and 6 J- in terms of SJ and SAG, we have 
equations that relate changes in the concentrations of 
the reactant and the product, c^ and Cp to change in the 
metabolic flux and free energy: 



SCr 


SJ+ 


SJ 




rSAG 










\ RT 


SCp 


SJ^ 


SJ 




/SAG 


Cp 






4 


\ RT 



(19) 



If the flux J changes as the metabolic network moves 
from one steady state to another, the concentrations Cr 
and Cp change according to Eq. 1191 The total relative 
change in reactant and product concentrations associated 
with the reaction can be calculated as. 



SCr 



J+ + J- 



J 



SJ_ 

J 



fSAG\ 
V RT j 



Jl + Jl fSAG\' 



J2 \ RT J 



(20) 



We deflne total sensitivity as 



'dim 



V = 



din J 



d In Cp 
din J 



j-Of + (1 - j+ey 



(21) 
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where 9 = i.^TT') steady-state change in the 

AG in response to a change in J. It can be thought of 
as the nonhnear biochemical resistance of the reaction. 
With given J+ and J_, rj reaches its minimum 77* = 

1/2 



when 6* is 6** — 



J++J- 



This is the least 



sensitive, maximal robustness condition for the reaction, 
irrespective of how the reaction is situated in a network. 
This result suggests that in maintaining certain reactions 
in metabolic pathways near equilibrium, i.e., J <C J+, 
biological systems may tend to minimize the sensitivities 
of concentrations of certain key species to perturbations 
in flux. This insight may be used as a lead for identifying 
regulatory sites in metabolic systems. 

This results allows us to associate chemical thermody- 
namics with robustness in a biochemical network. How- 
ever, as we have seen in the previous section, it is not 
possible to increase or decrease robustness by increasing 
or decreasing the enzyme activity for a given reaction 
r ^ p. By increasing enzyme activity, J-|_, J_, and J 
will all increase in the same proportion. Hence the min- 
imum value of rj is not affected. A futile cycle, however, 
is capable of regulating the minimal 77. 

Several special cases are important, 
(i) If the reaction is near equilibrium, J+, J_ ^ J, then 
we have the approximate relationship between AG and 
J 



AG = i?T In ( 1 - — 



Hence, 9 « j^i^ + j^)- Similarly we have 9 « j-(l-y-). 

Substitute these into Eq. (jJU) we have rj « ^^'^ ■ In 
this regime, the total sensitivity of the concentrations to 
the flux, 77, is simply proprotional to the driving force 
AG. It is clear that the sensitivity is related to how far 
the systems is away from equilibrium. 
(ii) If the concentration Cp is clamped, then the (5 J_ = 
in Eq. ^ and SAG = -^SJ. Hence 6* = ^, and 
V — 'j^- The ratio -j^ is known as irreversibility of the 
reaction; for an irreversible reaction it is unity, and it 
approaches zero as the reaction approaches equilibrium. 
This result was first obtained in 12]. 

(Hi) By a similar argument, if the Cr is clamped then 9 
— -J- , and 7] = This is our Eq. igj . 

(iv) If the reaction is a control point for the flux, then 
the crossover theorem for unbranched reaction pathway 
[28L I29I I dictates that for SJ > 0, one has Scr > and 
Scp < 0. This yields -j^ < 9 < j- . Substituting this into 
Eq. , one can easily show that 77 has upper and lower 
bounds 



J 



+ J1 



J 



(22) 



IV. COMPARISON OF TWO MECHANISMS 

FOR SHIFTING APPARENT EQUILIBRIA: 
MACROMOLECULAR BINDING AND FUTILE 
CYCLING 

Our analysis of sensitivity regulation would not be 
complete without discussing macromolecular binding as a 
mechanism for shifting apparent equilibria 30, 31] . For 
example, consider the case where a protein P binds p 
but not r in the reaction r ^ p with equilibrium con- 
stant Keq- In this case, although the presence of P does 
not change the equilibrium constant between r and p, it 
does change the equilibrium concentration ratio between 
total p {p plus pP) and r [sJ : 



„eq 



„eq 

_P_ 

„e9 



l-f 



K,g {1 + Kb^jS) , (23) 



in which Ki, is the association equilibrium constant be- 
tween P and p. Thus it is possible to reduce the sensi- 
tivity r]p in Eq. Q by introducing a binding protein for 
intermediary metabolite p (or increase rjr with a binding 
protein for r) that essentially serves as a buffer. Hence, 
there exist at least two distinct mechanisms of modu- 
lating robustness and sensitivity in biochemical systems. 
One is based on allosteric binding, and another is based 
on reversible chemical modification via futile cycle with 
energy expenditure. This issue was raised as early as 
1971 |32] after the discovery of enzyme regulation by re- 
versible chemical modification in terms of protein phos- 
phorylation and dephosphorylation via a futile cycle l3i| . 
Specifically, Fischer asks ^] (i) "Why have organisms 
found it advantageous to develop separate mechanisms 
to control the activity of enzymes, namely, by noncova- 
lent (allosteric) changes in structure mediated by appro- 
priate effectors (binding), and by covalent modifications 
(via futile cycle) of the proteins?" and (ii) "Why are 
these two mechanisms, i.e., noncovalent allosteric reg- 
ulation and covalent modification via phosphorylation, 
usually superimposed on one another even though the 
changes in conformation resulting in either activation or 
inhibition are essentially the same?" 

The essential difference between the binding mecha- 
nism and a futile cycle is that reversible binding does not 
expend energy. Yet there are important tradeoffs associ- 
ated with macromolecular binding. Specifically, buffering 
control by macromolecular binding requires the effectors, 
say a protein, to be present in an amount approximately 
equal to that of the metabolite. The acting principle 
here is that significant binding is necessary to achieve a 
significant shift in the effective equilibrium. The futile 
cycle approach requires only a relatively small amount 
of enzyme to catalyze the cyclic reactions. From a con- 
trol systems perspective, if feedback information is tied 
to the concentration of some protein, then the former ap- 
proach is ideal. On the other hand, if it is necessary to 
achieve signal amplification where only a small number 
of copies of a protein are activated, a phosphorylation- 
dephosphorylation futile cycle is appropriate for the task. 
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In addition, the costs of binding regulation and the futile 
cycles regulation are different. The former requires sig- 
nificant amount of biosynthesis of effectors in advance, 
while the latter requires only a small amount of enzymes 
for the hydrolysis reaction. The latter consumes energy 
during the regulation while the former pays in advance 
in the biosynthesis. In engineering terms, this is an is- 
sue of material cost versus energy utilization, an issue of 
overhead versus operational costs. 



V. REGULATORY SENSITIVITIES AND 
STATISTICAL THERMODYNAMICS 

The present work suggests an important relation 
between the regulatory sensitivity and robustness of 
metabolic systems in cell and the thermodynamics of 
biochemical reactions. This rather unexpected quantita- 
tive relationship deserves further investigation, especially 
from a systems biology perspective. In this section, we 
provide some initial discussions for the subject. 



A. Reversibility of biochemical reactions 



large concentration fluctuations. This is indeed the case. 
The sensitivities introduced in Eqs. Q and ||SJ) are inti- 
mately related to the steady-state concentration fluctu- 
ations in the intermediary metabolites [s^. To see this, 
we note that from Eqs. (0)) and © ?7i = J (cofci — J)~ 

and ?72 = J k2+k-i ~ ' respectively. 

It is shown in Appendix C (Eqs. 1471 and EHl and the 

discussion there) that the cofci and cofcifc2/(fc2 + fc-i) can 
be expressed as: 



and 



(cofci 



Cofclfc2 

fc2 + fc- 



(Aei(O)Aci(O) 



(Ac2(0)Ac2(t)) 



dt. 



(24) 



(25) 



The right-hand-sides of Eqs. (|24|l and H25|l are the con- 
centration fluctuations. Thus, we see that smaller the 
fluctuations, smaller the sensitivities. Elf et. al. have 
reached a similar conclusion based on a linear approx- 
imation of a nonlinear, stochastic biochemical kinetics 

m. 



All chemical reactions are reversible, although in some 
cases the backward rate may be negligibly small. Approx- 
imating such cases as irreversible is often acceptable in 
kinetic analysis, but is problematic for a thermodynamic 
analysis as illustrated in • In the biochemical liter- 
ature a reaction usually is operationally considered irre- 
versible if J+/J_ > 5 [23. When J+/J- « 1, i.e., J+, J_ 
3> J = J+ — J- , a reaction is considered to be near equi- 
librium and there is a linear relation between its flux and 
the chemical potential difference AG = i?Tln(J+/J_) 
~ RTJ/J'^. Thus the equilibrium forward and reverse 
flux J'^ = J'^ (divided by RT) is the conductance of a 
biochemical reaction operating near equilibrium ^3|- 



B. Sensitivity and concentration fluctuations 

Because of the thermal agitations and the stochastic 
nature of molecular reactions, the concentrations of a 
biochemical sgecies fluctuate, in a test tube reactions 
and in cells 39]. Rigorously, the concentrations (such 
as ci and C2) discussed so far represent the mean values 
of the concentrations of these species. If one is able to 
measure the fluctuating concentration of a species as a 
function of time, say Ci{t) for species 1, then one may 
calculate the magnitude of the concentration fluctuation 
by /„°°(Aci(0)Aci(t))dt where Aci(i) = ci{t) - (ci) is 
the deviation of ci(t) from its mean value. The notation 
(• • • ) denotes average of the quantity in the brackets, and 
(Aci (0) Aci(t)) is known as time-correlation function of 
fluctuating ci |3^ . 

One might imagine that a biochemical species with 
high sensitivity of concentration will display similarly 



C. Futile cycle and heat dissipation 

While the major function of futile cycles in signal 
transduction seems to be improving the performance of 
information processing against noise, the functional roles 
of futile cycles in metabolic systems potentially include 
improving sensitivity or robustness of metabolite concen- 
trations. In addition to affecting sensitivity, it is possible 
that futile cycles play an important role in generating 
heart and regulating temperature [42l |. Hence it is im- 
portant to analyze metabolic regulations with a systems 
perspective. In particular the intriguing suggestion [T^ 
that the futile cycles are important components in obe- 
sity, in weight loss, and even in the so-called Atkins' diet 
[ia. I43 . l45j . deserves further investigation. It is timely 
to rethink the issues of nutrition, thermogenesis [4^ , and 
futile cycles with a molecular as well as modern systems 
biology and metabolic engineering approach 0j ■ 
The present work provides a thermodynamic basis for 
studying futile cycles, which is likely to be essential in 
such studies. 



VI. SUMMARY 

Spending energy to gain control is not a foreign concept 
in engineering. Since this strategy is hallmark of control 
engineering, it should not come as surprise that biological 
cells use energy in controlling metabolism, transcription, 
and translation. Thus futile cycles which utilize biochem- 
ical energy are not necessarily "futile" ; they likely serve 
as mechanisms of biochemical regulation. 
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It has become increasingly clear from our recent work 
that futile cycles play a unique and essential role in cellu- 
lar regulation and signal transduction in the form of pro- 
tein phosphorylation dephosphorylation (and GTPase). 
While the phosphate group serves as a structural signal 
in enzyme activation, the phosphorylation reaction also 
provides a source of energy. The energy expenditure in 
fact increases the accuracy pjl, sensitivity speci- 
ficity @, and robustness [Slssl of the cellular informa- 
tion processing, overcoming cellular internal noises from 
thermal fluctuations, small copy numbers, and limited 
affinities. 

However to date, the role of futile cycles in metabolic 
regulation has been less quantitatively understood. In 
this work we have shown that metabolic futile cycles 
shift the effective equilibrium constants for biochemical 
reactions, modulating the sensitivity and robustness of 
intermediate concentrations to changes in flux. By shift- 
ing the effective equilibrium so that a reaction is moved 
away from equilibrium, the stoichiometric sensitivity is 
increased. Shifting it in the other direction reduces the 



sensitivity and enhances robustness. The direction of the 
shift, and hence a putative physiological role for a given 
futile cycle, can be diagnosed from the direction of op- 
eration of the futile cycle. When a futile cycle drives a 
reaction in the forward direction (clockwise in Figs. 2b 
and 4a), then the sensitivity of the concentration of in- 
termediate species 1 to the steady-state flux is enhanced. 
When a futile cycle drives the reaction in the direction 
opposite the net flux (counterclockwise in Figs. 2b and 
4a), then the sensitivity is reduced and the intermediate 
concentration is made robust to changes in the flux. 
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VIII. APPENDICES 
A. A: Sensitivity with different output flux controls 

Consider the reaction 

^ 1 2 ^ (26) 

With Co fixed, there are several ways to control the out- 
put flux J, either (i) directly controlling flux J, or (m) 
by controlling the rate constant fca in the downsteam re- 
action, 2 3, or (iii) by controlling the concentration of 
downstream species, C3 in 2^3, or the enzyme for the 

reaction. We show here that all these three cases yield 
identical expressions for the stoichiometric sensitivity. 
The kinetic equations for reaction system in H26I) are 

^ = fcico - + /s2)ci + A:_2C2, (27) 

dc2 f ^ 

—1 = fc2Ci - A:_2C2 - < k3C2 (28) 

'^^ [ + k^3C3 

where the last term in Eq. (|28|l represents the three cases 
above. 



10 



Consider the dynamic equation for C2 with the general 
expression, 

dc'] 

— = fcico - +A:2)ci + fc_2C2, (29) 

^ = fc2Ci - A:_2C2 - J(C2, C3, fcs, fc_3). (30) 

In steady-state, we have 

fciCo - (/c_i + k2)ci + fc_2C2 — 0, (31) 
fc2Ci - fc_2C2 - J(c2,C3,fc3,fc_3) = 0. (32) 

Therefore, we have 

- (fc-i + fe) (5ci) + fc_2 ((^cz) = 0, (33) 

k2{Sci) - k-2{Sc2) - SJ{c2,C3,k3,k-3) = 0. (34) 

Hence, no matter how J is changed, by ^3, by C3, or fc_3, 
as long as it is not by C2, we have 

Solving this pair of albegraic equations yields 

dci _ dc2 _ k2 + 

dJ fc_i' 9J /c_ifc_2 



(37) 



B. B: A simple example of how futile cycle drives 
concentration ratio away from equilibrium 

This example is motivated by the classic work on ki- 
netic proofreading 0. Let's consider a 3-state kinetic 
system with A, B and C shown in Fig. 5. The equi- 



librium constant ^12 



\BY- 



and reaction between 



(38) 



B and C is coupled to an energy source AGde = 
RT In /nilVi^'^rJ^ • We shall denote pseudo- first order rate 

[D]k°k3Ki2 ^ 

constants fc2 = k^[D\ and fc„2 = k°_2[E], and 7 = f^^^ 

^ (l + Ki2)[C] ^ fc_3 / fc-27 + fc3 ' 

[A] ([A] + [B]) fc3 V fc3 + A:-2 

where i?nn7 = —AGde is the amount of energy 
pumped input to the reaction. 7 = 1 for a closed sys- 
tem in equilibrium. For large 7, i.e., the kinetic cycle 
goes clockwise, it is possible to have 7fc_2 ^ ^3 ^ fc_2 
and, 

i " ^ (^) » 

the expected equilibrium ratio between B and Ai. On 
the other hand, for small 7, i.e., the energy pumping is 
counter-clockwise, it is possible to have 7/0-2 ^ ^3 ^ 
k-2, and 



[C]^k_3 fk3\ fc-3 

[A] ^ ks \k.2J /C3 ■ 



(40) 



C. C: Concentration fluctuations in open systems 

Here we consider concentration fluctuations in the ki- 
netic pathway with constant source cq and sink C3 = 0: 



ki k2 i-o 

4 1 4 2 — 

fc-l fe-2 



3. 



(41) 



For ageneral theory on open, linear biochemical networks 
see [25|. In stochastic terms, the probability of the num- 
bers of species 1 and 2 being m and n at time t, P{m, n, t) 
satisfies the chemical master equation |^ |^ 



dp{m, n, t) 
dt 



(fciCo -I- TO(fc_i + k2) + (fc_2 + k3)n)p{m,n) 



+kiCop{m—l, n) + (m+l)fc_ip(m+l, n)+{m+l)k2p{m+l, n—1) 



+ {n + l)fc_2P(™ - 1, n + 1) + (n + l)fc3p(m, n + 1) (42) 

From Eq. 142(1 it is easy to show that the mean values of 
m and n follow the standard deterministic kinetic equa- 
tions 



dt 
d 

dt 



(m) = fciCo - (A:_i -|- k2){m) + fc_2(n) (43) 



{n) = k2{m) - (fc_2 + k3){n) 



(44) 



{{Amy) = fciCo + (fc_i +/C2)(™) +fc_2(") 

-2(fc_i + k2){{A'mf) + 2fc_2(AmAn), 



and furthermore their variances and covariance, 

d 
di 

^((Anf) = k2{m) + ik^2 + k3) (n) 
dt 

+2k2{AmAn) - 2(fc_2 + k3){{Anf), 
— (AmAn) ^ -k2{m) - k^2{n) + k2{{Amf) + k^2{{Anf) 
-(/c_i + k2 + fc_2 + fc3)(AmAn), 

where Am = m — (m)^^ and An = n — (n)^^. The steady 
state mean values 



fcifc2Co 



fc-l(fc-2 + ^3) + ^2^:3 ' fc_i(fc_2 + ^3) + ^^3 ' 

(45) 

which agree with Eqs. ^ and Q if we substitute J — 
(ri)^*fc3, and more, 

{{A7n)y = {my% {{An)y = (r^)''^ (AmAn)"" 0. 

(46) 

In fact, one can verify that the stationary probability dis- 
tribution for m and n is Poissonian [2^. ^^"^ 

''^"Ll ^^"^ ^-'^^ stationary distribution for Eq. H42|l . 
However, while Eq. H4t)|l indicates that the steady-state 
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fluctuations for m and n are independent. It does not 
mean that m and n are compleltely independent. In fact, 
we have (Am(O)An(t)) 7^ for t 7^ 0: 

(Am(O)ATO(i)) = 

(m) 



(-(A2 + + k2)e^^' + (Ai + fc_i + fc2)e^=*) 

(47) 



Ai — A2 

(An(O)An(t)} = 



{nY 



Ai - A, 



(-(A2 + fc-2 + h)e^'' + (Ai + fc-2 + fc3)e^^*) 

(48) 

(49) 



(ATO(O)An(i)) - (e/i* 



(An(O)Am(i)} 



Ai — A2 



-,A2t\ 



Ai — A2 

where Ai and A2 are the two eigenvalues of the matrix 



(50) 



-fc_l - /C2 
k2 



k-2 
-k-2 - h 



(51) 



Note that the concentation ci (02) and number of 
molecules m (n) differ by the volume of the system. 
Hence, dentifying <aM£#) ^ W£W*)> 



and i£if£#l 



= ''"(nj"i)*i^ ; ^-iid integrating Eqs. ITTIandBSIwe arrive at 
Eqs.E|andE3 

We shall also point out that the stochastic dynamics on 
a mesoscopic scale can be described by a Fokker-Planck 
equation |50ll51| 



d 



P(rn, n, t) 



dt 

in which 

r(m, n) 



(52) 



A 
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FIG. 5: A realistic mechanism, motivated by the kinetic 
proofreading theory for protein biosynthesis @ , showing how 
a driven cyclic reaction, i.e., futile cycle, regulates the popula- 
tion ratio between A and C away from their equilibrium ratio 
k^a/kg. The chemical driving force for the cycle is —AGde = 
RT\n{Ki2kik3[D]/{kUk-^3\E])). 



representing two species and three reactions in Eq. 

According to the Fokker-Planck equation, the sta- 
tionary distribution P^^im^n) is Gaussian centered ar- 
round {(m)'''^ , ("■)''*) with covariant matrix cr satisfying 
Ik, 31, 531 



F(T + aF^ = T {{mY% {nY') ■ 



(55) 



It can be easily verified via matrix multiplication that the 
solution to the Eq. (|55|l is precisely Eq. H46|l . We should 
point out that the fluctuation covariance is not FF^ 
as recently suggested '54]. The information on the con- 
centration fluctuations is not contained in the relaxation 
rate {F) alone; it has to depend on T which represents 
the rate of fluctuations p^ . 



(fciCo + (fc-l + k2)m + fc_27^- — fc2TO — fc_2^ 
— fc2m — k^2n k2m + (fc_2 + fcs)'^ 

(53) 

and F is given in Eq. H51() . We see that T can be ex- 
pressed as ^S{J'^ + J^)S'^ where the stoichiometric ma- 
trix [Hi 



